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摘 要: 结合 博 斯 腾 湖 1960 一 2018 年 水 位 .出 入 湖 径 流 以 及 气象 站 点 实测 资料 ,采用 集合 经 验 模 态 分 解 (Ensem- 
ble Empirical Mode Decomposition ,EEMD ) .水 量 平衡 和 气候 弹性 方法 ,对 近 60 a 博 斯 腾 湖 水 位 变化 及 其 影响 因素 进 
行 了 详细 分 析 。 结 果 表 明 :(1) 1960 一 2018 年 博 斯 腾 湖 水 位 总 体 呈 下 降 态 势 ,具体 表现 为 "下 降 - 上 升 - 下 降 - 上 升 ” 
四 个 阶段 。(2) 在 年 际 尺 度 上 水 位 存在 准 3~4a、 准 8~9 a 的 周期 性 振荡 ,而 年 代 际 尺度 上 表现 出 准 29~30 a 和 准 33~ 
34 a 的 周期 性 变化 ,(3) 1960—2018 年 降水 气温 和 潜在 蒸 散 发 对 开 都 河 、 黄 水 沟 和 下 并 径流 的 累积 贡献 率 分 别 达 
85.1% .42.1% Fil 23.8% ,而 下 执 面 、 其 他 气象 变量 和 人 为 等 因素 累积 对 径流 的 贡献 率 分 别 约 为 14.9%、57.9% 和 
76.2%。(4) 对 不 同 阶段 博 斯 腾 湖 水 位 变化 原因 分 析 :1960 一 1987 年 水 位 急剧 下 降 的 主要 原因 同人 湖 径流 减少 和 湖 
面 蒸 发 量 大 有 关 ; 和 气温 升 高 和 降水 量 增加 导致 人 湖水 量 增加 是 1988—2002 年 水 位 显著 升 高 的 主要 原因 ;和信 湖 径流 


减少 和 出 湖水 量 增多 ,导致 2003 一 2014 年 水 位 显著 下 降 ; 博 


是 2015 一 2018 年 水 位 明显 上 升 的 主要 原因 。 


斯 腾 湖 入 湖水 量 的 显著 增加 及 对 出 湖水 量 的 严格 控制 


关键 词 : 博 斯 腾 湖 ; 水 位 变化 ; 集合 经 验 模 态 分 解 ; 水 量 平衡 ;气候 弹性 ; 影响 因素 


干旱 区 内 陆 湖泊 是 对 全 球 气候 变化 响应 最 为 
敏感 的 地 理 单元 之 一 , 易 受 人 类 活动 的 影响 一 。 湖 
泊 水 位 过 低 和 过 高 都 会 之 来 不 利 影响 ,其 中 ,水 位 
下 降 过 低 将 导致 湿地 面积 萎缩 ,植被 退化 ,生物 多 
样 性 减少 和 渔业 损失 。 而 水 位 过 高 ,水 面 面 积 扩 
大 ,将 加 剧 土 壤 盐 碱 化 和 农业 损失 ,增加 粮食 风险 ， 
而 超过 一 定 生态 水 位 还 会 产生 一 系列 洪水 灾害 , 造 
成 一 定 的 社会 经 济 损失 ““"。 湖 油水 位 变化 的 原因 
非常 复杂 ,尤其 是 干旱 地 区 的 平原 湖泊 , 易 受 其 地 
理 特征 、 气 候 变 化 、 冰 川 覆 盖 和 人 类 活动 等 因素 的 
影响 ”。 因 此 ,加 强 干 旱 区 湖泊 变化 的 科学 研究 对 
深入 了 解 其 水 量 平衡 具有 重要 意义 ,并 为 如 何 调节 
这 些 重要 水 资源 提供 线索 。 

干旱 区 内 陆 湖 变 化 最 易 受 气候 变化 和 人 类 活 
动 的 影响 ,其 中 ,平原 区 尾 间 湖 由 于 气候 变化 和 
人 为 影响 而 急剧 收缩 ,高 山 湖泊 由 于 降水 变化 较 小 
而 保持 相对 稳定 , 否 吐 型 湖泊 由 于 河流 入 湖 口 用 水 
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量 的 变化 , 比 其 他 两 种 类 型 的 湖泊 变化 更 灵活 *。 
黄 坤 等 ”采用 交叉 小 波 研究 了 各 种 气象 因素 对 巴 
尔 喀什 湖水 位 变化 的 影响 。Liu 等 Jing SE" Al Li 
等 ”利用 遥感 数据 研究 新 疆 不 同 区 域 湖泊 退缩 和 
扩张 ,认为 湖泊 变化 是 自然 和 人 为 因素 共同 作用 的 
结果 。 苏 向 明 等 “采用 高 分 遥感 影像 分 析 艾 比 湖 
的 年 内 和 年 际 变 化 特征 ,认为 年 内 湖泊 变化 主要 影 
响 因 素 为 农业 灌溉 用 水 。 阿 依 努 尔 * 买 买 提 等 “ 采 
用 统计 分 析 方 法 对 伊 塞 克 湖水 位 变化 的 影响 因素 
进行 探讨 分 析 , 认 为 气候 变化 为 其 水 位 变化 的 第 一 
影响 因素 ,人 为 因素 为 第 二 影响 因素 。 作 为 国内 最 
大 的 内 陆 淡 水 湖 一 一 博 斯 腾 湖 ,其 水 位 变化 一 直 是 
许多 学 者 关注 的 热点 问题 Yao 等 和 Guo 
等 “采用 水 量 平衡 和 气候 弹性 方法 ,结果 显示 不 同 
时 期 博 斯 腾 湖 水 位 变化 的 影响 因素 明显 不 同 。 朱 
长 明 等 “通过 采用 多 源 和 遥感 数据 研究 博 斯 腾 湖 水 
位 面积 和 水 量 时 间 动 态 变 化 过 程 与 演变 特征 , 认 
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为 其 水 资源 呈现 波动 减少 主要 受气 候 变 化 和 人 类 
活动 的 共同 影响 。 本 文 基于 FEEMD 方法 ”0 和 水 量 平 
衡 原 理 ,详尽 分 析 了 博 斯 腾 湖 1960 一 2018 年 水 位 变 
化 趋势 以 及 在 不 同时 期 的 非 线性 变化 过 程 , 采 用 和气 
候 弹性 方法 ,探讨 了 博 斯 腾 湖 水 位 在 不 同 阶段 对 水 
文 气象 因子 的 响应 程度 ,分 析 了 湖泊 水 位 对 气候 变 
化 和 人 类 活动 的 响应 过 程 , 旨 在 为 博 斯 腾 湖 流域 水 
资源 可 持续 利用 和 生态 环境 可 持续 发 展 提供 科技 
支撑 。 


1 研究 区 概况 与 研究 方法 


1.1 研究 区 概况 

博 斯 腾 湖 (41°56’'N~42°14'N,86°40'E~87°56'E) 
地 处 新 疆 南部 ,是 中 国 最 大 的 内 陆 淡 水 湖 (图 1)。 
该 地 区 是 典型 的 内 陆 气候 , 博 斯 腾 湖 湖区 多 年 平均 
气温 约 6.3 % ,年 均 降 水 量 仅 约 70 mm ,潜在 蒸发 量 
高 达 2000 mm25。 博 斯 腾 湖 作为 新 疆 典 型 的 知 吐 型 
湖泊 ,既是 塔里木 河 四 源流 之 一 一 一 开 都 河 的 尾 闻 
湖 ,同时 也 是 下 游 孔 汰 河 的 补给 湖 。 博 斯 腾 湖 水 源 
来 自 天 山南 坡 的 开 都 河 、 黄 水 沟 .清水 河 及 诸 小 河 
流 , 其 中 开 都 河和 黄 水 沟 作 为 博 斯 腾 湖 的 主要 补给 
源流 ,分别 占 总 人 湖水 量 的 84.6% 和 6.3%5 。 博 斯 
腾 湖 作为 开 都 一 孔 八 河流 域 水 资源 时 空 分 布 的 天 
然 调节 库 , 同 时 也 是 新 疆 巴 音 郭 楞 蒙古 自治 州 地 区 
(简称 巴 州 ) 的 重要 水 资源 储存 库 ,具有 调节 河川 径 


84°E 


流 净化 水 质 等 生态 功能 ,被 誉 为 巴 州 人 民 的 “母亲 
湖 ”。 
1.2 数据 资料 的 选取 

(1) 径流 数据 。 大 山口 和 宝 浪 苏 木 ( 开 都 河 )、 
黄 水 沟 ( 黄 水 沟 ) 和 塔 什 店 ( 孔 省 河 ) 水 文 站 的 
1960—2018 实测 年 径流 量 数据 由 新 疆 塔里木 河流 
域 管理 局 提供 。 

(2) 气象 数据 。 巴 音 布鲁克 、 巴 伦 台 和 看 普 气 
象 站 1960—2018 年 的 日 气温 降水 风速、 蒸发 H 
照 时 数 .水汽 压 数据 来 自 中 国 气象 数据 服务 网 (http: 
//data.cma.cn/) o 

(3) 湖水 位 数据 。 博 斯 腾 湖 1960—2018 年 实 
测 年 水 位 数据 和 1986—2018 年 实测 月 水 位 数据 来 
源 于 博 斯 腾 湖 管理 局 。 

(4) 土地 利用 数据 。1980、1995、2005 年 和 
2015 年 土地 利用 数据 来 源 于 中 国 科学 院 资源 环境 
科学 与 数据 中 心 (http://www.resdc.cn/data.aspx?DA- 
TAID=197) ,分 辩 率 为 1 km 的 栅 格 数据 ,土地 利用 
类 型 包括 6 个 一 级 类 型 以 及 25 个 二 级 类 型 。 
1980 、1995 、2005 年 各 期 数据 的 遥感 解 译 主要 使 用 
了 Landsat-TM/ETM 遥感 影像 数据 ,而 2015 年 土地 
利用 /覆盖 数据 更 新 主要 使 用 Landsat 8 遥感 影像 
BE 
13 研究 方法 
13.1 EEMD 方法 ”集合 经 验 模 态 分 解 (EEMD ) 
PRON ,适用 于 非 线性 和 非 平 稳 序 列 信号 分 析 ,能 够 


86°E 88°E 
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图 1 博 斯 腾 湖 流域 示意 图 
Fig. 1 Location of the Bosten Lake Basin 
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干旱 区 研究 
很 好 地 提取 趋势 和 周期 信息 ”。 在 竺 分析 的 原始 AQ=AQ,+ AQ, (7) 


言 号 序列 中 二 加 上 给 定 振幅 的 白 噪 声 序列 ,将 加 入 
白 噪 声 后 的 信号 进行 EMD 分 解 ( 分 解 为 本 征 函 数 
IMFs) ,将 多 次 分 解 得 到 的 IMFs 进行 集合 平均 ,使 加 
入 的 白 噪 声 互 相抵 消 , 并 将 其 作为 最 终 的 分 解 结果 : 
区 的 = t500 (1) 
式 中 : CO 为 原始 信号 经 过 EEMD 分 解 变 换 后 得 到 
的 第 7 个 IMF 分 量 ;NN 为 白 噪声 增加 数 ; CO) 为 第 i 
次 加 入 白 噪 声 后 分 解 所 得 的 第 /个 IMF 分 量 。 
1.3.2 M-K ttk Mann-Kendall 检验 2 的 统计 


tt SH: 
5=¥ $ sma- 2) 
统计 UF, 的 序列 值 为 : 
aD pa (3) 
JVar(S,) 

检验 统计 区 域 的 均值 和 方差 如 下 : 
ECS )= 1) (4) 
Vary DY Gs 


通过 对 开 都 河 BEL ACTA) AES E 8 AEE FT 8 
检验 ( 表 1) ,本 人 研究 将 时 间 序 列 分 为 突变 前 基准 期 
(1960 一 1996 年 ,1960 一 1992 年 ,1960 一 1991 年 ) 和 
突变 后 (1997 一 2018 4 , 1993—2017 4F , 1992—2016 
年 ) 两 个 时 期 。 

1.3.3 气候 弹性 方法 根据 Zheng 等 ”提出 的 气候 
弹性 方法 ,计算 开 都 河和 黄 水 沟 径 流 对 气候 变化 的 
敏感 性 : 

ZTC-9) 


=a = (6) 
'Q JOA- 


R1 博 斯 腾 湖 水 位 和 不 同 水 文 站 的 径流 
Mann-Kendall 趋 势 及 突变 检验 
Tab.1 Mann-Kendall trend and mutation test of 
Bosten Lake water level and runoff from different 


hydrological stations 
因子 Z Sig ”a=0.01 突变 时 间 趋势 ”显著 性 
水 位 -4.24 0.00002 2.56 19664 PRE 显著 
大 山口 径流 3.19 0.00140 2.56 19964F EF 显著 
黄 水 沟 径流 2.86 0.00430 2.56 199248 “上升 “显著 
Aue 0.48 0.62990 2.56 199148 ”上升 不 显著 


注 : Z 表示 Z 分 数 ,表示 变化 的 显著 性 ; Sig 表示 双边 检验 的 P 值 。 


AQ, =(epAP/P+ es AET JET, +erAT/T)O (8) 
式 中 : 6 为 气象 因子 的 弹性 系数 (i 表示 气温 7 , 降 
水 P ARTERE ET, ); O, 和子 分 别 是 气象 因子 
(如 降水 气温 等 ) 引 起 的 年 径流 变化 量 和 不 同 气 候 
变量 的 值 , O 是 径流 量 的 变化 比例 , 2 是 平均 径流 ， 
为 气象 因子 平均 值 ; AQ 、AP AET, 和 A7 分 别 
为 流域 内 径流 、 降 水 、 潜 在 菩 散 发 和 气温 的 变化 量 ; 
AO.、AO, 分 别 是 气象 因子 和 其 他 因素 (包括 下 热 
面 特性 和 人 为 因素 等 ) 引 起 的 径流 变化 。 
1.3.4 水 量 平 衡 ” 博 斯 腾 湖 是 吞吐 型 湖泊 ,其 水 量 
平衡 方程 可 以 表示 为 二 
AV=(O,,+P,- On, -E + AV, (9) 
式 中 : AV 是 博 斯 腾 湖 水 量变 化 量 ; On , Oo, 是 入 出 
湖 的 水 量 ; P, 为 湖面 降雨 量 ; 5, 湖面 蒸发 量 , AV, 
为 误差 和 地 下 水 交换 量 引 起 的 湖水 变化 量 。 
1.3.5 潜在 蒸 散发 的 计算 ” 洪 在 蒸发 数据 是 根据 气 
象 站 的 逐日 气象 观测 资料 ,采用 联合 国 粮农 组 织 推 
荐 的 Penman-Monteith 公式 计算 所 得 ”| 
0.408A(R, - G)+»y—200 
ET= 


T2737) 
A +y(1 + 0.34u,) 

式 中 : ET NAMA; R, 为 净 辐 射 ; G 为 土壤 热 通 
Ht; e, 是 饱和 蒸汽 压 ; e, 是 实际 水 汽 压 ; 了 是 地 面 
以 上 2 高度 处 空气 温度 ; u, 是 地 面 以 上 2 mm 高 处 
的 风速 ; A 为 饱和 水 汽 压 与 温度 曲线 上 的 斜率 ; y 是 
湿度 计算 常数 。 这 些 参数 是 由 平均 日 最 高 和 最 低 
温度 ,平均 温度 ,实际 水 汽 压 ,平均 风速 ,日 照 时 数 ， 
平均 相对 湿度 数据 和 其 他 经 验 指 标 直 接 计算 得 到 。 
1.3.6 湖面 蒸发 量 计算 “” 博 斯 腾 湖 湖面 的 年 蒸发 量 
采用 Zhou 等 计算 方法 ,可 以 表示 为 : 

E,=Ey Xkxs (11) 
式 中 : 妃 是 湖面 单位 面积 的 年 蒸发 量 ; Fe 是 口径 
为 20 mm 的 蒸发 严 测 量 的 蒸发 量 ; 大 是 湖面 蒸发 量 
修正 系数 ; s 是 湖面 面积 。 


2 结果 分 析 


2.1 博 斯 腾 湖 水 位 变化 

根据 Mann-Kendall 趋势 检验 结果 , 近 60 a 博 斯 
腾 湖 水 位 整体 呈 显 著 下 降 趋 势 (P<0.01, 表 1) ,具体 
表现 为 四 个 阶段 (图 2a):1960 一 1987 年 以 0.08 m-a” 


(10) 
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1980 2000 2018 
年 份 


1046.70 


(b) 年 内 
1046.65 | 


1046.60 


1046.55 上 


水 位 /m 


1046.50 


1046.45 


1046.40 
12345 67 89101112 


月 份 


图 2 1960 一 2018 年 博 斯 腾 湖 水 位 年 际 变 化 (a) 和 年 内 各 月 水 位 变化 (b) 
Fig. 2 Inter-annual (a) and monthly (b) changes of water level in Bosten Lake form 1960 to 2018 


的 速率 下 降 了 近 3 m,1988 一 2002 年 以 0.26 ma 的 
速率 上 升 了 近 4 m,2003 一 2014 年 水 位 以 0.18 mea! 
的 速率 下 降 了 近 3.2 m, 自 2015 年 以 来 , 博 斯 腾 湖 水 
位 以 0.72 m a :速率 上 升 近 1.9m( 图 2a)。 其 中 ， 
1987 年 和 2002 年 分 别 是 博 斯 腾 湖水 位 最 低 (1045 
m) 和 最 高 (1049.36 m) 的 两 个 年 份 ,相差 4.36 mo 

博 斯 腾 湖 水 位 年 内 变化 呈现 双 峰 型 形态 (图 
2b)。 其 中 9 月 的 水 位 最 高 (1046.65 m) ,1 月 最 低 
(1046.45 m)。4 一 9 月 是 一 年 中 高 水 位 时 期 ,11 月 
至 型 年 1 月 是 一 年 中 低 水 位 时 期 。 春 季 气 温 回 升 ， 
促使 上 游 山 区 冰川 积 雪 消 融 加 速 , 进 而 使 得 下 游 径 
流量 增加 ,湖泊 水 位 抬升 ,至 4 月 水 位 达 1046.59 
m。4 月 以 后 ,由 于 博 斯 腾 湖 水 量 调度 和 周边 农业 活 
动 的 影响 ,湖泊 水 位 有 所 下 降 。 夏 秋季 降水 量 增 
加 ,季节 性 冰川 积 雪 完全 消融 ,入 湖 径 流量 增加 , 故 
9 月 份 湖泊 水 位 最 高 。 冬 季 开 都 河上 游 来 水 量 少 ， 
年 内 径流 量 最 小 ,湖泊 处 于 低 水 位 状态 。 
2.2 水 位 的 多 尺度 变化 特征 

图 3 是 对 博 斯 腾 湖 1960—2018 年 水 位 距 平 进 
行 EEMD 分 解 得 到 4 个 IMF 分 量 和 1 个 趋势 分 量 
(RES); FEP ,分解 所 得 到 的 下 Ms 依次 反映 了 原始 
水 位 距 平 从 高 频 到 低频 不 同时 间 尺 度 的 波动 变化 ， 
趋势 项 RES 表示 1960 一 2018 年 水 位 随时 间 变 化 的 
整体 演变 趋势 。 结 合 表 2,IMF1 未 通过 5% 的 显著 性 
检验 ,表示 准 3~4 a 周期 属于 弱 周 期 ;IMF2、IMF3 和 
IMF4 均 通过 了 5% 的 显著 性 检验 ,说 明 IMF2、IMF3 
和 IMF4 分 量 非常 显著 ,水 位 具有 准 8~9 a HE 29~30 
a 和 准 33~34 a 的 主 周期 性 振荡 。 

各 分 量 方差 贡献 率 显示 , 博 斯 腾 湖 水 位 年 际 变 
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图 3 博 斯 腾 湖 1960 一 2018 年 水 位 距 平 各 IMF 分 量 及 趋势 项 


Fig. 3 The IMF and trend components of water level 
anomaly for Bosten Lake from 1960 to 2018 


表 2 博 斯 腾 湖 水 位 距 平 各 IMEF 分 量 的 周期 、 方 差 贡 献 
率 及 其 与 径流 距 平 序列 的 相关 性 
Tab. 2 Periods and variance contribution rates of IMF 
components and their correlations with the water level 


anomaly series for Bosten Lake 


IMF 分 量 IMF1 IMF2 IMF3 IMF4 RES 
周期 /a 3.28 8.89 29.50 33.83 
与 水 位 距 平 的 相关 系数 0.29 0.23 0.90" 0.39" 0.29" 
贡献 率 /% 6.88 7.61 83.04 1.68 0.79 
显著 性 检验 /% <90 >95 >95 >95 


注 :*x# * 分 别 在 0.01 .0.05 水 平 上 显著 相关 。 


化 约 占 14.49% , 而 年 代 际 变化 为 85.51% ,年 代 际 变 
化 对 水 位 整体 变化 的 贡献 高 于 年 际 变化 。 年 代 际 
振荡 在 水 位 长 期 变化 中 占 主 导 ,结合 水 位 各 分 量 
(MFs) 与 原始 水 位 距 平 的 相关 关系 ,发 现年 代 际 变 
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化 同 原始 径流 距 平 呈 非常 显著 的 正 相 关 。 结 合 图 3 
和 表 2,IMF3 表征 的 准 29~30a 周 期 性 振荡 的 方差 贡 
献 率 最 大 , 约 为 83.04% ,该 周期 尺度 在 1980 年 之 后 
有 明显 振荡 ;IMF1 表示 的 是 准 3~4 a 周期 振荡 ,主要 
JE JE 2000—2005 年 和 2010 年 之 后 表现 的 比较 明 
显 ,其 对 水 位 距 平 总 体 方差 贡献 率 为 6.88%;IMF2 
代表 的 准 8~9 a 的 周期 性 振荡 ,方差 贡献 率 为 
7.61% , 其 周期 性 振荡 在 整个 研究 时 期 均 有 表现 ; 
IMF4 表示 的 是 水 位 准 33~34 a 的 周期 变化 ,其 方差 
贡献 率 最 小 (1.68% ) ,在 此 时 间 尺 度 上 , 博 斯 腾 湖 水 
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位 在 1980 一 1990 年 振幅 处 于 较 低 水 平 ;由 趋势 项 变 
化 可 以 看 出 , 博 斯 腾 湖水 位 在 1960 一 2018 年 间 总 体 
上 表现 出 非 线 性 下 降 趋 热 。 
2.3 博 斯 腾 湖 水 量 平衡 

研究 发 现 , 博 斯 腾 湖 湖面 蒸发 量 和 出 湖水 量 对 
水 位 下 降 的 贡献 率 分 别 达 到 了 47% 和 58% ,由 图 4 
和 表 3 可知 ,1960 一 1987 年 和 人 湖水 量 和 湖面 降水 量 
均 呈 下 降 趋势 ,而 该 时 期 博 斯 腾 湖 水 位 下 降 了 3.04 
m, 湖 水 量 减少 了 约 3.10 km’;1988 一 2002 年 , 博 斯 腾 
湖人 湖水 量 和 湖面 降水 量 对 水 位 的 变化 贡献 率 分 
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4 博 斯 腾 湖 水 量 平 衡 要 素 湖 面 蒸发 量 (a) 湖面 降水 量 (b) 和信 湖 径 流量 (c) 和 出 湖 径 流量 (d) 的 变化 趋势 
Fig.4 Changes of water balance elements in Bosten Lake: Evaporation of lake surface (a), precipitation on lake surface (b), runoff 
into the lake (c) and discharge from the lake (d) 


3 博 斯 腾 湖 四 个 时 期 的 水 量 平衡 


Tab.3 Lake water balance in Bosten Lake during four periods 


aa 湖泊 补给 量 湖水 消耗 量 地 下 水 交换 
湖泊 水 量变 化 一 一 一 一 一 a = == TZ te A 
阶段 BCom at) ， 人潮 径流 量 ke ”湖面 降水 量 占 比 WERE 占 比 ”出 湖水 量 be 量 + 误差 占 比 /% 
(km-a)  /% Ikmsat)  /% /aa /% Ikma) /% /km a ) 
1960—1987 4E -0.11 1.92 48 0.07 2 0.93 23 1.15 29 -0.07 2 
1988—2002 年 0.28 2.72 55 0.10 2 0.83 18 1.48 30 -0.23 5 
2003—2014 4E -0.28 2.14 46 0.07 1 0.94 20 1.65 35 -0.11 2 
2015—2018 4E 0.50 2.84 53 0.08 1 1.08 20 1.45 27 -0.09 2 
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别 为 115% 和 4% ,而 该 时 期 入 湖水 量 .湖面 降水 量 、 
湖面 蒸发 量 和 出 湖水 量 均 呈 上 升 趋势 ,同期 湖水 位 
上 升 了 约 4.01 m, 水 量 增加 4.13 km°s 2003—2014 
年 ,湖面 蒸发 量 和 出 湖水 量 对 湖泊 水 位 下 降 贡 献 量 
分 别 为 43% 和 75% ,研究 发 现 该 时 期 入 湖水 量 ,湖面 
降水 量 ,湖面 蒸发 量 和 出 湖水 量 均 呈 下 降 趋 势 , 同 
期 湖水 位 下 降 约 3.24m ,水量 减少 了 约 3.41 km SR 
Ti , 2015—2018 年 , 博 斯 腾 湖 人 湖水 量 增加 ,湖面 降 
水 量 下 降 ,湖面 蒸发 量 和 出 湖水 量 增加 ,而 同期 入 
湖水 量 和 湖面 降水 量 对 水 位 的 变化 贡献 率 分 别 达 
112% 和 3% ,使 得 同期 水 位 上 升 1.93 m 左 右 , 水 量 增 
加 了 约 1.98 km’, 

2.4 博 斯 腾 湖 水 位 变化 的 影响 因素 分 析 

2.4.1 气象 因子 变化 对 径流 的 影响 “径流 是 博 斯 腾 
湖 补给 的 主导 因素 ,从 图 5 可 以 看 出 , 博 斯 腾 湖 水 位 
变化 同 开 都 河和 黄 水 沟 径 流 变化 高 度 一 致 (P< 
0.05) , 故 其 径流 变化 将 直接 影响 博 斯 腾 湖 水 位 ,而 


气温 .降水 和 潜在 蒸 散发 主要 通过 影响 其 上 游 径流 
进而 对 下 游 湖水 位 产生 影响 。 降 水 .气温 和 潜在 薰 
散发 与 径流 的 相关 性 系数 ppo、pro 和 Pero 表明 
( 表 4) , 开 都 河和 黄 水 沟 的 降水 与 径流 的 相关 性 最 
显著 ,表明 降水 是 影响 径流 变化 的 最 关键 因子 。 根 
据 表 4 中 降水 (P) 气温 (7) 潜在 蒸 散发 (ETu) 和 径 
流 (O) 的 阶段 变化 ,利用 气候 弹性 方法 计算 得 到 开 
都 河源 流 区 的 气候 弹性 系数 ep = 0.65, e, = -0.31 和 
sir=-0.67( 表 5), 即 1996 年 之 后 开 都 河流 域 降水 、 
气温 和 潜在 薰 散 发 每 增加 /上 升 10% ,其 对 应 的 径流 
量 将 增加 6.5% .3% 和 -6.7%。 黄 水 沟 流域 径流 的 气 
候 弹 性 系数 sp= 0.93, ep = 0.98 和 sor=-2.85, 即 3 
个 气象 要 素 变 化 10% , 将 会 引起 径流 分 别 发 生 
9.3% 、9.8% 和 -2.85% 的 变化 。 丰 着 气象 站 和 水 文 站 
位 于 焉 着 盆地 腹地 ,水文 站 以 上 的 集 水 区 面积 约 为 
4134.2 km ,在 该 区 域 , 干 旱 指 数 ( p=E7wWP ) 大 于 


1.2 T T 
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图 5 标准 化 后 的 博 斯 腾 湖 水 位 变化 与 主要 人 湖 河 流 径流 变化 


Fig. 5 Changes of Bosten Lake water level and runoff of two rivers into the lake after standardization 


4 不 同时 期 气象 因子 和 径流 的 变化 


Tab.4 Changes of climatic elements and runoff in different periods 


阶段 Pimm ET/mm TIC Q/mm R, o Pro Pro Pero 

山口 1960 一 1996 年 257.30 632.60 -4.70 174.70 0.68 2.46 0.58” 0.32 -2.44 
1997—2018 年 303.90 642.90 -3.60 211.50 0.70 2.12 0.61” 0.05 0.52" 

1960—2018 4 274.70 636.50 -4.30 188.40 0.69 2.32 0.69" 0.43" -0.17 
黄 水 沟 1960—1992 年 192.52 1048.32 6.17 57.66 0.30 5.45 0.60" 0.21 -0.47" 
1993—2017 4 237.57 977.91 7.14 83.56 0.35 4.12 0.80" -0.01 -0.42° 
1960—2017 年 211.94 1017.97 6.59 68.83 0.32 4.80 0.75" 0.34 -0.57” 

KE 1960—1991 4E 71.23 919.20 8.20 55.80 0.78 12.90 0.18 0.20 -0.08 
1992—2016 4Ẹ 82.95 864.00 9.10 63.80 0.77 10.40 0.17 0.13 -0.54" 
1960 一 2016 年 76.38 895.01 8.60 59.33 0.78 11.72 0.21 0.22" -0.36" 


注 :#*#* # 分 别 在 0.01.0.05 水 平 上 显著 相关 。 
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10 ,径流 系数 (R.=OP) 小 于 0.8( 表 4) ,属于 典型 的 
干旱 区 水 文 气候 特征 ,其 径流 的 气候 弹性 系数 为 
sp=0.13、sr=0.90 和 szr=-0.98, 即 当 降 水 量 增加 
10% 时 ,其 相应 的 年 径流 量 将 增加 1.3% ; 而 气温 升 
高 10% 时 ,将 导致 年 径流 量 增加 9% ; 当 洪 在 蒸 散发 
减少 10% 时 ,年 径流 量 将 增加 9.8%。 

由 表 5 可 知 ,三 个 水 文 站 的 降水 弹性 系数 ep = 
0.65 .0.93 和 0.13 ,相对 于 基准 期 ,降水 分 别 增 加 了 
46.61 mm .45.05 mm 和 11.72 mm ,根据 公式 (8) 可 以 
计算 出 大 山口 1997—2018 年 BE 7K YA 1993—2017 
EVI R Fa E 1992—2016 年 径流 深 比 基准 期 分 别 增 
加 约 20.56 mm、12.30 mm 和 1.20 m; 气 温 弹 性 系数 
e, 分 别 为 -0.31、0.67 和 0.9, 温 度 分 别 上 升 1.10 C, 
0.97% 和 0.90% ,引起 大 山口 RAIA ALES Be 
深 相 对 于 基准 期 分 别 增加 12.67 mm、9.00 mm 和 
3.10 mm; 潜在 蒸 散 发 弹性 系数 © pr, 分 别 为 -0.67、 
-2.85 和 -0.98 ,三 个 集 水 区 的 潜在 蒸 散发 量 分 别 增 
加 了 10.33 mm -70.41 mm 和 -55.20 mm, 导致 径流 
深 相 对 于 基准 期 分 别 减 小 约 1.91 mm 10.40 mm 和 
6.20 mm, 

研究 发 现 ,降水 .气温 和 潜在 蒸 散发 的 综合 影 
响 下 ,1997 一 2018 年 开 都 河 年 经 流 深 增加 了 31.32 
mm, 占 径流 深 增加 量 的 85.1% ,而 其 他 因素 作用 下 
导致 年 经 流 深 增加 5.50 mm, 占 流量 变化 的 14.9% 。 
1993 一 2017 年 黄 水 沟 气 象 因子 变化 累积 引起 径流 
深 增加 10.90 mm , 约 占 径流 增加 总 量 的 42.1% ,而 其 


他 因素 变化 导致 径流 深 增加 15.00 mm , 约 占 径流 变 
化 量 的 57.9%。 而 在 农业 比较 发 达 的 在 弟 舍 地 人 类 
活动 比较 频繁 ,1992 一 2016 年 气温 .降水 和 潜在 蒸 
散发 变化 直接 引起 的 径流 深 总 的 变化 量 约 为 1.90 
mm, 约 占 径 流 总 变化 量 的 23.8% , 远 小 于 其 他 因素 
所 导致 的 径流 深度 变化 (76.3% )( 表 5)。 这 说 明 , 相 
对 于 基准 期 ,1990s 之 后 ,入 湖 径流 除了 对 气温 、 降 
水 和 潜在 蔡 散 发 变化 敏感 外 ,其 他 要 素 对 径流 变化 
也 很 重要 。 
2.4.2 人 为 因素 对 水 位 变化 的 影响 “研究 表明 , 宝 
浪 苏 木 人 湖水 量 分 别 占 1960 一 1987 年 .1988 一 2002 
年 .2003 一 2014 年 和 2015 一 2018 年 大 山口 径流 量 的 
60% .70% 、59% 和 72% , 耗 水 量 分 别 约 占 大 山口 径 
流 总 量 的 36% .26% .27% 和 17%( 表 6) ,这 也 意味 着 
流域 内 的 引 耗 水 量 对 博 斯 腾 湖 水 量变 化 有 很 大 的 
影响 。 

总 的 来 看 ,1988 一 2002 年 博 斯 腾 湖 流域 开始 进 
入 丰 水 期 2002 年 博 斯 腾 湖 流域 大 山口 径流 量 达 到 
近 60a 的 最 大 值 (4.49 km’) , 博 斯 腾 湖水 位 也 达到 最 
高 。 根 据 已 有 资料 后,1950s 流域 内 农田 灌溉 面积 
约 为 100 km? ,至 1960 年 农田 灌溉 面积 已 经 增长 至 
578 km ,大 量 的 耕地 开垦 以 及 水 资源 的 肆意 利用 ， 
导致 1987 年 之 前 入 湖 径 流量 不 断 减 少 。1986 年 之 
后 ,由 于 灌溉 定额 制度 的 实施 ,水 资源 的 利用 效率 得 
以 提高 ,流域 内 灌溉 用 水 量 得 到 有 效 的 控制 ”?” 。 
1980 .1995 .2005 .2015 年 博 斯 腾 湖 流域 土地 利用 变 
化 数据 显示 ( 表 7) ,流域 内 耕地 面积 仍然 在 不 断 增 


表 5 气温 、 降 水 和 潜在 蒸 散 发 变化 对 径流 的 影响 


Tab.S Effects of changes in temperature precipitation and potential evapotranspiration on runoff 


弹性 系数 AP AET, AT AQ APIP AET, jy ATIT AQIQ AQP) AQT) AQ(ET) AQ. Ad, 

Ep A Em /mm /mm /mm /mm 1% ET, 1% 1% /mm /mm /mm /mm /mm 
大 山 0.65 -0.31 -0.67 46.61 10.33 1.10 36.82 18.11 1.63 -23.40 21.08 20.56 12.67 -1.91 31.32 5.50 
黄 水 沟 0.93 0.98 -2.85 45.05 -70.41 0.97 25.90 23.00 -7.00 16.00 45.00 12.30 9.00 -10.40 10.90 15.00 
Be 0.13 0.90 -0.98 11.72 -55.20 0.90 8.00 16.00 -6.00 11.00 14.00 1.20 3.10 —6.20 1.90 6.10 


表 6 大 山口 径流 量 、 宝 浪 苏 木 入 湖 径流 量 与 大 山口 至 宝 浪 苏 木 段 引 耗 水 量变 化 
Tab. 6 Changes of runoff in Dashankou and runoff of Baolangsumu into Bosten Lake and water consumption from 


Dashankou to Baolangsumu Station 


时 间 阶 段 大 山口 径流 量 /km 宝 浪 苏 木 人 湖 径流 量 /km 大 山口 一 宝 浪 苏 木 段 引 耗 水 量 /km 
1960 一 1987 年 3.22 1.92 1.16 
1988—2002 年 3.86 2.72 1.02 
2003—2014 4 3.64 2.14 0.97 
2015—2018 年 3.94 2.82 0.66 
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表 7 博 斯 腾 湖 流域 土地 利用 变化 


Tab.7 Land use changes in Bosten Lake basin 


年 份 耕地 /km 林地 /km 草地 /km 水 域 /km 未 利用 地 /km 城乡 居民 工矿 用 地 /km 
1980 1691 428 20143 1133 19325 78 
1995 1775 450 19796 1134 19536 100 
2005 1900 495 19619 1162 19486 128 
2015 2207 455 19636 1162 19202 133 


加 ,至 2015 年 达到 约 2207 km’, 1980—2015 年 期 间 ， 
博 斯 腾 湖 流域 耕地 .城乡 居民 用 地 面积 均 有 所 增 
加 ,其 中 耕地 增长 率 高 达 约 30.5%。 耕 地 和 生产 生 
活用 地 的 不 断 增加 ,会 导致 流域 内 下 垫 面 发 生变 
化 ,可 能 会 对 流域 内 气候 和 径流 产生 重要 影响 。 
2015 年 之 后 ,水 位 快速 回升 的 主要 原因 除了 入 湖 径 
流量 的 增加 外 ,参考 博 斯 腾 湖 适宜 生态 水 位 ,人 
为 对 出 湖水 量 的 严格 控制 也 是 主要 因素 。 

由 于 水 资源 的 无 限制 利用 ,塔里木 河 下 游人 遭受 
了 严重 的 生态 退化 ,为 了 恢复 塔里木 河 下 游 的 “ 绿 
EER” ,政府 于 2000 年 实施 了 生态 输 水 工程 ,将 水 
从 博 斯 腾 湖 输送 到 大 西海 子 水 库 , 其 所 用 的 水 主要 
来 自 博 斯 腾 湖 和 节 水 灌溉 ,这 进一步 加 剧 了 博 斯 腾 
湖 流域 水 资源 短缺 。 从 2002 一 2010 年 , 博 斯 腾 湖 的 
总 生态 输 水 量 为 1.91 km 中 。 这 也 是 造成 2003 一 
2014 年 间 博 斯 腾 湖 水 位 下 降 的 原因 。 


3 讨论 


开 都 河 径 流 对 气象 因子 的 敏感 性 分 析 结 果 表 
明 , 降 水、 气温 和 潜在 蔡 散 发 对 径流 的 累计 贡献 率 
约 为 85.1%, 其 他 因素 (包括 人 为 和 气象 因子 间 相 互 
作用 导致 径流 的 变化 ) 对 径流 的 贡献 率 约 为 14.9%， 
这 同 Chen 等 3 利用 水 文敏 感性 方法 分 析 得 到 开 都 
河 径流 对 气候 变化 和 人 类 活动 的 响应 结果 基本 一 
致 ,Li 等 中 在 阿克苏 河流 域 采用 气候 弹性 方法 的 研 
究 也 得 出 类 似 的 结论 。 降 水 气温 和 潜在 蒸发 对 径流 
的 贡献 率 差 别 较 大 ,降水 明显 要 高 于 气温 和 潜在 蒸发 
的 影响 ,Wang 等 ”在 漆 河 流域 也 得 出 类 似 的 结论 。 

1960—1987 年 .1988 一 2002 年 2003—2014 年 
和 2015 一 2018 年 博 斯 腾 湖 水 位 变化 的 主要 原因 分 
别 为 :径流 量 减少 和 湖面 蒸发 量 大 气候 变化 引起 
入 湖 径 流量 增加 、 入 湖 径 流量 减少 和 出 湖水 量 增 
多 、 人 为 对 出 湖水 量 的 控制 和 入 湖 径 流量 的 增加 ， 
这 表明 径流 是 引起 水 位 变化 周期 的 最 主要 因素 。 
通过 对 气象 因子 的 EEMD 分 析 结 果 可 知 ,在 近 60 a 


时 间 尺 度 上 年 际 和 年 代 际 尺度 上 的 降水 与 径流 呈 
显著 正 相 关 关 系 ,而 气温 和 潜在 蒸 散发 与 径流 主要 
在 年 代 际 尺度 上 显著 相关 ,表明 降水 是 影响 径流 变 
化 的 最 主要 的 气象 因素 。 本 文采 用 了 小 波 分 析 方 
法 进一步 判断 EEMD 的 分 解 效果 ,结果 显示 28 a 周 
期 对 小 波 方差 的 贡献 率 最 大 ,与 EEMD 分 解 结果 在 
年 代 际 尺度 上 基本 一 致 ,但 与 郭 梦 京 等 ”采用 连续 
小 波 分 析 求 得 水 位 变化 周期 结果 有 所 不 同 , 可 能 是 
由 于 研究 的 时 间 段 不 同 所 致 ,但 相 比较 而 言 ,本 人 研 
究 的 数据 时 间 序 列 更 长 ,结果 的 可 靠 性 更 高 。 开 都 
河流 域 大 山口 径流 在 年 际 尺 度 上 存在 准 3~4 a HE 
7~8 a 的 周期 性 振荡 ,在 年 代 际 尺度 上 存在 准 13~14 
a 和 准 29~30a 的 周期 性 变化 。 降 水 存在 准 3~4a、 准 
7~8 a 的 年 际 尺度 变化 , 准 14~15 a YE 51~52 a 的 年 
代 际 变化 。 湖 面 蒸发 量 在 年 际 尺度 上 呈 准 2~3 a HE 
9~10 a 的 周期 变化 ,在 年 代 际 尺度 上 有 准 29~30 a, 
准 39~40 a 的 周期 性 。 宏 观 上 ,这 些 不 同 尺度 的 周 
期 还 可 能 来 自 于 大 气 环流 以 及 太阳 活动 的 准 周 期 
变化 5 , 受 西风 环流 及 北冰洋 气流 影响 的 博 斯 腾 湖 
流域 ”, 水 位 变化 可 能 存在 一 定 的 关联 。Wang 等” 
分 析 区 域 气候 和 遥 相 关 对 博 湖 流域 水 文 过 程 的 影 
响 表明 , 博 斯 腾 湖 流域 水 文 过 程 受 北半球 极 区 指数 
涡 (VPA) 和 西藏 高 压 (XZH ) 的 影响 (统计 上 ) 显 著 ， 
尤其 是 气温 和 径流 。VPA 的 变化 导致 博 斯 腾 湖 流 
域 中 纬度 西风 酸 和 山脊 系统 发 生变 化 ,进而 影响 博 
斯 腾 湖 流域 的 气候 。 

根据 入 湖 量 与 湖面 蒸腾 蒸发 量 的 比值 , 博 湖 的 
蒜 发 蒸腾 损耗 对 于 水 位 变化 有 着 很 重要 的 影响 ( 表 
3)。 湖 面 的 芋 发 量 随 着 湖面 积 增 大 而 增加 ,尤其 是 
在 气温 较 高 的 月 份 。 因 此 ,可 以 通过 协调 湖泊 水 位 
和 湖面 蒸发 量 之 间 的 关系 ,减少 水 资源 损失 ,提高 
水 资源 利用 率 ”。 已 有 人 研究 表明 '» 开 都 河 年 径流 
量 主要 集中 在 夏季 (6 一 8 月 ) , 约 占 全 年 径流 量 的 
44% ,气温 升 高 加 速 了 冰川 积 雪 的 融化 ,降雨 频 率 
较 高 ,径流 大 幅 增 加 ,发 生 洪涝 灾害 风险 的 概率 增 
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大 '”。 因 此 ,在 调 水 时 间 分 配 上 ,可 以 在 主要 农 作 
WEKE 湖水 蒸发 最 强 的 5 一 9 月 向 下 游 孔 淮河 流 
域 调 水 ,同时 ,每 年 10 月 至 翌年 4 月 作物 收获 后 , 温 
度 和 蒸发 速率 较 低 , 博 斯 腾 湖 水 位 在 冬季 和 春季 保 
持 较 高 的 水 位 ,从 而 最 大 限度 地 减少 博 斯 腾 湖 的 年 
蒸发 量 损失 ,提高 水 资源 利用 率 2”。 

除 气候 变化 外 ,人 类 活动 对 博 斯 腾 湖 水 位 变化 
也 有 一 定 的 影响 ,如 农业 灌溉 引水 ,水 利 工程 建设 
等 (中 。 自 1958 年 以 来 ,人 类 活动 虽然 仍 小 于 气候 变 
化 对 湖水 位 变化 带 来 的 影响 ,但 对 博 斯 腾 湖 水 位 变 
化 的 影响 呈 增 加 趋势 ”1, 主要 体现 在 人 类 活动 频 
繁 的 出 山口 以 下 平原 和 盆地 的 主要 入 湖 径 流 消 耗 
区 的 影响 “'。 博 斯 腾 湖 流域 水 资源 系统 是 脆弱 的 ， 
气候 变化 加 剧 了 不 稳定 水 资源 的 不 确定 性 ,但 人 为 
用 水 因素 对 博 斯 腾 湖 的 未 来 同样 起 着 关键 作用 。 


4 4 it 

(1) 1960 一 2018 年 博 斯 腾 湖 水 位 总 体 呈 下 降 态 
势 ,经 历 了 下 降 - 上 升 -下降 - 上 升 的 过 程 。1960 一 
1987 年 以 约 0.08 m+ a! HJ IE F BE, 1988—2002 年 
以 0.26 m- a! 的 速率 上 升 ,2003 一 2014 年 水 位 以 
0.18 ma 的 速率 下 降 ,2015 年 以 来 ,水 位 以 0.72 me 
a' 速 率 上 升 , 近 60a 间 , 博 斯 腾 湖 水 位 总 共 下 降 了 
0.45 mo 

(2) 博 斯 腾 湖 1960 一 2018 年 水 位 的 EEMD 分 解 
结果 表明 ,水 位 呈 准 8~9 a、 准 29~30 a 和 准 33~34 a 
的 周期 性 振荡 。 近 60 a 来 水 位 变化 在 年 际 尺度 上 存 
在 准 3~4 aE 8~9 a 的 周期 性 振荡 ,在 年 代 际 尺度 上 
存在 准 29~30 a 和 准 33~34 a 的 周期 性 变化 ,其 中 在 
29~30 a 时 间 尺 度 上 周期 振荡 比较 明显 。 

(3) 气象 因子 对 径流 敏感 性 分 析 表 明 ,1997 一 
2018 年 开 都 河流 域 期 间 ,降水 .气温 和 潜在 蒸 散发 
对 径流 贡献 率 分 别 为 57%、34% 和 5%。1993 一 2017 
年 黄 水 沟 流 域 降水 .气温 和 潜在 车 散发 对 径流 量 的 
贡献 率 分 别 为 47% .35% 和 40%。1992 一 2016 4E SG 
普 盆 地 气温 降水 和 洪 在 蒸 散发 变化 直接 引起 的 径 
流 变化 量 约 占 径流 总 变化 量 的 15% .39% Fl 78% © 

(4) 近 60 a 不 同时 期 博 斯 腾 湖 水 位 变化 来 看 ， 
1960—1987 年 水 位 下 降 主要 同和 湖 径 流 减 少 
(-0.11 km-a”) ,湖面 蒸发 耗 水 量 大 (1.02 km +a! A 
关 ; 1988—2002 年 降水 增加 ,气温 升 高 和 湖面 潜在 
蒸 散发 量 较 小 ,导致 径流 量 增 加 (0.28 km? +a!) , 引 耗 


水 量 下 降 (-0.006 km*+a') ,是 湖泊 水 位 上 升 的 主导 
素 ;2003 一 2014 年 比 1988 一 2002 年 平均 人 湖水 量 
减少 约 0.58 km ,出 湖水 量 增加 0.11 km( 生 态 输 水 
工程 为 主 ) ,是 该 阶段 湖水 位 下 降 的 主要 原因 。 
2015 年 以 来 ,入 湖水 量 增加 以 及 人 为 对 出 湖水 量 的 
合理 控制 ,水 位 持续 上 升 至 1047.59 m(2018 年 )。 
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Analysis of the change in water level and its influencing factors 
on Bosten Lake from 1960 to 2018 
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Abstract: Using measurements of the water level and runoff into the lake as well as data from meteorological 
stations of Bosten Lake from 1960 to 2018, this paper provides a detailed analysis of the trend in water level 
change based on ensemble empirical mode decompsition and the water balance model. In addition, this paper 
uses the climate elasticity method to explore the response of the water level to hydroclimatic factors in different 
time periods, and then analyzes the complex response of the lake water level to climate change and human 
activities. The results revealed the following: (1) The lake level shows a significant decreasing trend from 1960 to 
2018, which is shown in the four periods as “descending-rising-descending-rising”. (2) Ensemble empirical mode 
decomposition analysis results show that the water level of Bosten Lake has periodic oscillations of 3-4 a and 8- 
9 a on the interannual scale, and 29-30 a and 33-34 a on the interdecadal scale. (3) The contributions of precipitation, 
temperature, and potential evapotranspiration to runoff to the Kaidu River, Huangshui River, and Yangqi Basin 
from 1960 to 2018 were approximately 85.1%, 42.1%, and 23.8%, respectively. Other factors contributed 
approximately 14.9%, 57.9%, and 76.2% to the runoff. (4) The prime reason for the decline in water level in 
Bosten Lake from 1960 to 1987 was the decrease in runoff into the lake and substantial evaporation of the lake 
surface. Increase in lake water intake due to the climate change was the main cause of rise in water level from 
1988 to 2002. The decrease in runoff into the lake and increase in the amount of lake water are responsible factors 
for the decreased water level from 2003 to 2014. An increase in lake water intake and strict controls on discharge 
from 2015 to 2018 were the major reasons for the rising water level in Bosten Lake. 

Keywords: Bosten Lake; water level change; ensemble empirical mode decomposition; water balance; climate 
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